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Abstract. A new approach is proposed to the analysis of generalized synchronization 
of multidimensional chaotic systems. The approach is based on the symbolic analysis 
of discrete sequences in the basis of a finite T-alphabet. In fact, the symbols of the 
T-alphabet encode the shape (the geometric structure) of a trajectory of a dynamical 
system. Investigation of symbolic sequences allows one to diagnose various regimes 
of chaos synchronization, including generalized synchronization. The characteristics 
introduced allow one to detect and study the restructuring and intermittency behavior 
of attractors in systems (the time structure of synchronization). The measure of T- 
synchronization proposed is generalized without restrictions to complex ensembles of 
strongly nonstationary and nonidentical large-dimensional oscillators with arbitrary 
configuration and network (lattice) topology. The main features of the method are 
illustrated by an example. 
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1. Introduction 


Synchronization is one of the fundamental concepts of the theory of nonlinear dy¬ 
namics and chaos theory. This phenomenon is widespread in nature, science, engineering, 
and society [1]. One of important manifestations of this phenomenon is the synchronization 
of chaotic oscillations, which was experimentally observed in various physical applications 
(see Si and references therein) such as radio oscillators, mechanical systems, lasers, elec¬ 
trochemical oscillators, plasma and gas discharge, and quantum systems. The study of this 
phenomenon is also very important from the viewpoint of its application to information 
transmission [^, cryptographic coding with the use of deterministic chaotic oscillations, 
and quantum computation [^|^. 

There are several types of synchronization of chaotic oscillation s p| : generalized syn- 

lag synchroniza- 


chronization 
tion 
tion fl5 


9 , complete synchronization 10 , antisynchronization 


frequency synchronization phase synchronization 14 , time scale synchroniza- 
and T-synchronization pTl. For each type, an appropriate analytic apparatus and 
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diagnostic methods have been developed. Nevertheless, intensive investigations are being 
continned that are aimed, on the one hand, at the examination of different types of syn¬ 
chronization from nniqne positions and, on the other hand, at the search for new types of 
synchronons behavior that do not fall nnder the above-mentioned types. In spite of the long 
history of the stndy of synchronization of chaotic oscillations, many important problems in 
this held remain unsolved. 

These include generalized synchronization in the form 


y(f) = F[x(t), r]. 


( 1 . 1 ) 


where x and y are multidimensional synchronized systems, F is a function of generalized link 
between the systems, and r is a delay vector between the phase variables of the systems x 
and y. 


In this paper, we develop an original method for the diagnostics and quantitative 
measurement of the characteristics of generalized synchronization of chaotic systems, which 
is aimed at the integrated study of the time structure of synchronism through the analysis 
of the so-called T-synchronization 


16,17 


The method is based on the formalism of symbolic CTQ-analysis proposed by the 
18,^ (the abbreviation CTQ stands for three alphabets with which the 
C, T, and Q). 


present author 


method operates: C, T, and Q). One should note that symbolic dynamics, for all of its 
seeming external simplicity, is a very strongly substantiated tool for the analysis of nonlinear 
dynamical systems 


20 22 


This article is an expanded version of the report 23 


2. The Symbolic CTQ-analysis 


Denote a discrete dynamical system as a mapping 


Sfc+i f (sfc, p) 


( 2 . 2 ) 


with the following properties: s G S C M”, /cgKCZ, pGPC n G 1, A^, m G 1, M. 
In formula (2.2), s is a state variable of the system and p is a vector of parameters. With 


mapping (2.2), we associate its trajectory in the space S x K, which has the form of a 
semisequence fc G 1, iP. 


2.1. T-alphabet 


Dehne the initial mapping, which encodes (in terms of the hnal T-alphabet) the shape 
of the n-th component of the sequence [^,19 : 


„(D „(D „(D \ . — \T°“^\ rpaip\ rpaipi 1 /r) q\ 

Sfc-n Sfc \n, -[Ik \l,---,Ik \n,---,Ik IvJ • 


The graphic diagrams illustrating the geometry of the symbols for the k-th 

sample and the n-th phase variable are shown in Figure [Tj 
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Figure 1. Geometry of T-alphabet sym¬ 
bols. 



Figure 2. Transition matrix of T- 
alphabet symbols. 


Strictly speaking, the mapping (2.3) is defined by the relations: 
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and 

As+ = 

I 

= 

(n) 

t-l-1 

(n) 







(2.4) 


Here As_ = 

Thus, the T-alphabet includes the following set of symbols: 

= iTo, Tl, T2, T3N, T3P, T4N, T4P, T5N, T5P, 

T6S, T6, T6L, T7S, T7, T7L, T8N, T8P}. (2.5) 


One can see from (2.5) that the symbol T^'^\n is encoded as Ti, where i is the right- 


hand side of the symbol codes of the alphabet T“‘^. In turn, the symbol is encoded in 
terms of Tii ■ ■ ■ in ■ ■ ■ iN, see (2.3). The full alphabet T“‘^|iV, which encodes the shape of 


the trajectory of the multidimensional sequence consists of 17^ symbols. 
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2.2. Q-alphabet 


In addition to the symbols we introduce the symbols 

n“‘^l = _v I D“‘^l D"'^l 1 

Vfc In — k \n -tfc+ilnj Vfc [Vfc II5 • • • 5 Vfc Inj • • • 5 Vfc |ArJ 


( 2 , 6 ) 


-.Qipi 


All admissible transitions constitute a set of symbols of the alphabet 3 |„. 

These transitions are shown in Figure 

The symbol Q^'^ln is encoded as Qij, where i and j are the right-hand sides of 
the symbol codes of the alphabet T"*^ for the states k and A; -|- 1, respectively. In turn, 
the symbol is encoded in terms of QA ■ ■ ■ in ■ ■ - ji ■ ■ ■ jn ■ ■ ■ Jn, see (2.6). The full 

alphabet which encodes the shape of the trajectory of the sequence consists 

of 107^ symbols (see Figure]^. 


2.3. Symbolic TQ-image of a dynamical system 

Let us introduce a directed graph 

YTQ ^ ^ yr ^ x 7 |Ar, C Q 7 |iV, 


(2.7) 


which is a complete symbolic TQ-image of the dynamical system (2.2). By dehnition, is 


the vertex set and E'" is the edge set of F^'^. According to its topology, the graph F^*^ has 
no multiple arcs but has loops. 

By analogy with ( 2.7| ), introduce a directed graph 

= (V^U, E^n) , C T7, E^n C Q7, (2.8) 

which is a particular symbolic TQ-image of the dynamical system with respect to its n-th 
phase variable. Denote the graph T'^^\n corresponding to the full alphabets T""^ and Q""^ 
by FJQ. 


A particular symbolic TQ-image (2.8) can be obtained from the graph F^*^ by the 


gluing (or identihcation) of its vertices, redirection of the edges incident to these vertices, 
and the removal of the arising multiple edges. 

Introduce the operation of gluing: 


Let n e V^, e e V^, v' e V^l^, and e' e V^|„. Then 

F = Ti, V = lii ■ ■ ■ in ■ ■ ■ in, 

e = Qi j, c = Qii ■ ■ ■ in ''' iN ji ■ ■ ■ jn ■ ■ ■ Jn- 

The gluing of vertices and edges is formally expressed by the following conditions: 


3 Tj : TA 


iN e V^, in = j, 


E |n 3 Q i J . Q ki • • • kn ■ ■ ■ kN h ■ ■ ■ ■ ■ ■ In ^ E , kn ii In T 


(2,9) 


( 2 , 10 ) 


( 2 . 11 ) 


The graph (2.8) can be weighted (on its vertices and edges) by the occurrence fre¬ 
quency of characters * in the sequence 


A* I = 

^ \n 


|M* 


[JM* 


0 ^ < 1 , 


( 2 . 12 ) 
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where | • | is the cardinality of the set and * is a symbol of which the multiset M*|„ consists: 


3 : T“^|n\T = *, * G T7\T, (2.13a) 

: M*U 3 Q7\n : Qrin\Q = *, * G Qr\Q. (2.13b) 


Similar characteristics A"*" and A*^ are determined for the graph Note that the calcu¬ 
lation of A"*" and A*^ allows one to quantitatively assess various properties of the trajectory 
of the sequence {s/d^yin the space S x K, including the Markov characteristic of the 


sequence {T7}7i 20 


22 


3. T-synchronization 


Let us make a remark. For simplicity, but without loss of generality, suppose that 
the time sequence of dimension N is formed by a combination of the phase variables 

of N one-dimensional dynamical systems; i.e., suppose that is the value of the phase 
variable of the nth system at the fcth instant of time. 


3.1. Complete T-synchronization 

Definition. Dynamical systems are synchronous at time instant k in the sense of 


Complete T-synchronization 17 if the condition = 1 is satished, where 


Jh. 


1 Tn, = ...=T“>’i„=...=Tru, 

0 otherwise. 


(3.14) 


Thus, {is the indicator sequence of T-synchronization. 


Taking into account possible antisynchronization 11 between the systems, we should 
also consider all possible variants of inversion of their phase variables: —>■ — 1 ■ 8^**^ In 

this case, for the nth system, a change of symbols in the fcth sample occurs according 

to the scheme 

TO <H- TO, 


T1 GG T2, T3N GG T5P, T3P GG T5N, T4N GG T4P, 
T6S O T7S, T6 O T7, T6L O T7L, T8N O T8P. 


(3.15) 


Denote the variants of inversion by number m. The total number of variants of inversion 
is M = 2^-F 


Synchronization between systems can also be set in the lag regime 12 . To detect this 


synchronization, one should move a little the phase trajectories of the systems with respect 
to each other (the shifts are ^ 0): 


(\ . rjiCH-^p I rjiCtip \ . rjiOiip I rjiOL^^p \ . rjiOiip I "1 

|l -^fc+hi hi • ■ ■ 1 \n ^ J-k+hjn, ■ ■ ■ , J-k \N ^ J-k+hNlN f ■ 


(3.16) 


The antisynchronization and lag synchronization regimes may coexist; therefore, when 
calculating a partial integral coefficient of synchronism, we take this fact into consideration: 


K* 


^ ^ 5^ M {m, h}. 


k=k* 


(3.17) 
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where k* = 1 + max (hi, ..., hiv), K* = K + min (hi, ..., Hn), and K is the length of the 
sequence 

On the basis of the partial coefficient, we calculate the total integral coefficient of 
synchronism of the systems: 


= maxmaxh^hJ 0 ^ ^ 1, (3.18) 

m h ’ 

i.e., we take a combination of shifts between the trajectories of the systems and a variant 
of inversion of their phase variables that, taken together, provide the maximum number of 
samples k satisfying the condition J*, = 1. 


4. Time structure of synchronization of chaotic systems 


The quantity h* introduced in (3.18) characterizes the synchronism of the systems 
on average over a period of tx — ti. As mentioned in the Introduction, most investigations 
on the synchronization of chaos are usually restricted to this situation. However, often a 
researcher may be interested in the time structure of synchronization of systems. Recall that 
by this structure one means the spikes in the synchronous behavior of the phase variables of 
the systems between which the synchronism level is characterized by a small quantity, i.e., 
one means intermittent behavior . 


In 17 , the present author introduced the concept of a synchronous domain SD 


a set of samples of a time series that satisfy the condition (V is the symbol of the logical 
operation OR) 


SDr : {t/fc' = 1, Jk” = 0 y k" = 0, Jkiii = 0 V k'" = A" + 1} , 
k' e foSD, foSD + Lr-l, k" = - 1, k'" = bf^ + 


(4.19) 


where 6^^, and r are the emergence time, the length, and the ordinal number of 

a synchronous domain, respectively. In this case, the following conditions are satished: 

^ K, and the total number of synchronous domains (in the original sequence) ^ 
(AT l)div2. 

To quantitatively describe the structure of synchronization of systems, the author 
introduced in 17 the spectral density function of synchronous domains SD: 




i]. iel, A-, 


(4.20) 


r=l 


where (5[-, ■] is the Kronecker delta. 

To analyze the degree of degeneracy of the structure of synchronous domains, we 
additionally define a quantity ~ the entropy of the structure of synchronous domains 
(according to Shannon) which makes sense for 5^ > 0: 

pSD [p] ^ [L] ^4 21) 

*=1 E [^] 

i=l 

It follows from Shannons entropy properties that the entropy is minimal (A®° = 0) 
when the spectrum 77®^ [L] is degenerate (all synchronous domains have the same length) and 
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Figure 3. Basic characteristics of the time structure of synchronism. 


maximal in the case of a uniform comb spectrum H^'^[L] with the maximum 

number of different lengths of synchronization domains equal to 


= min 


yJl + S5^K - 1 


, K-5^K + l\ , E^^ = \xiW^ 


cmb'! 


(4.22) 


where [aj is the integer part of a. 


On the basis of (4.21) and (4.22), we dehne the relative entropy of the structure of 
synchronous domains: 

^SD 

(4.23) 


A|° = 


^SD 


It makes sense to apply the quantity A|P when the researcher should compare synchroniza¬ 
tion cases with different values of 5^ and/or K. 

Nevertheless, for the full description of the intermittent behavior of chaotic systems 
during synchronization, it is obviously insufficient to study only synchronous domains SD. 
To obtain a complete and closed idea of the time structure of synchronism of dynamical 


systems (a complete and closed representation of the intermittency structure), in 16 the 
present author introduced the concept of a desynchronous domain SD - a set of samples of 
a time series satisfying the condition 


SD^ : {t/fc' — 0, Tfc" — 1 V k" — 0, Jk'" — 1 V k"' — K -|- I} , 
k' e + Lf^ _ k'' = hf^-l, k''' = hf^ + Lf, 


(4.24) 


where bf.^, T®, and r are the emergence time, the length, and the ordinal number of a 
desynchronous domain SD, respectively. 

The meaning of the characteristics introduced in this section is demonstrated in Fig¬ 
ure 1^ (see for additional information). 


4.1. Generalized T-synchronization 


mzation 


11 


It follows from the dehnition of the synchronization condition (3.14) that the ana¬ 
lyzer proposed evaluates the complete synchronization level 10 and detects antisynchro- 
with lag synchronization (I^ precisely in the alphabetic representation T"'^. 


However, according to the dehnition of the geometry of the symbols of the T-alphabet (2.4), 


complete synchronization at the level of the samples of is a wider phenomenon com¬ 
pared with the complete synchronization at the level of - the samples of the sequence 
itself. The T-synchronism of dynamical systems (with respect to the set of phase variables s) 
is considered from the viewpoint of the shape (geometric structure) of the trajectories of the 
systems in the extended phase space. By the shape (geometric structure) of a trajectory of a 
dynamical system in the extended phase space is meant its certain invariant under uniform 
translations and dilations of the trajectory in the space of phase variables. 
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Thus, in a sense, the T-synchronization deals with the topological aspects of synchro- 

Hence, this opens a possibility for the application of 


nization of dynamical systems 20,21 


the analyzer proposed to the study of generalized synchronization of chaos [9]. 

To this end, we introduce two additions that relax the requirements imposed in Sec¬ 
tion 3.1 on the complete T-synchronization. The hrst is the rejection of the equality of 


symbols in (3.14), and the second is the rejection of the maximization of 5* in (3.18). 


The rejection of the equality of symbols in (3.14) allows one to proceed to the following 
dehnition. 

Dynamical systems are synchronous at time instant k in the sense of Generalized 
T-synchronization if the condition = 1 is satisfied, where 


Jh 


1 e M-^T, 

0 otherwise. 


MJT C T7|iV. 


(4.25) 


The structure of the set is not trivial. First, the cardinality of the set is 

|mJT| =niin||vni|, IV^nl, 

Second, the condition 


\N\ 


(4.26) 


MJ^liJ^l, nel,N, 


(4.27) 


is always satisfied. 

As a rule, conditions (4.26) and (4.27) correspond to different variants of constructing 
the set Let 3 be the set of all admissible variants of constructing the 

set where j is the number of a variant. 

The number of variants of constructing the set is 


Pw-l N-l 


2=0 n=l 




\N\- 


Note that the indices n in (4.28) are rearranged as follows: 


Ivnil ^ ^ Iv^J ^ ... ^ Iv^^l 


If the condition 


IV'^llU ...= |V^Ll= ... = |V^|Ar| 


(4.28) 


(4.29) 


is valid, then 




Pjv —1 

n(^’»- 


i=0 


N-l 



T^{Pn) 
r(i + Piv)’ 


where T is the gamma function. 

It should also be noted that, in the general case, represents an upper estimate, 
because the matrix of A"*" may contain zero elements; i.e., not all possible combinations of 
elementary symbols are realized. 

Thus, complete T-synchronization is a special case of generalized T-synchronization. 
Figure illustrates the structure of the set for A = 2 in two cases of complete T- 
synchronization, the direct synchronization and antisynchronization. 
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Figure 4. The structure of the set M'^'^ for iV = 2 in two cases of complete T-synchronization: direct 
synchronization (green) and antisynchronization (orange). 


According to (4.28), in the case of generalized T-synchronization, the problem arises 
of choosing a set from among the family that is optimal with respect to some 

criterion. In the general case, this is a combinatorial optimization problem. Let ns introduce 
a generalized algorithm 


^ max, G 

where is a objective function. 


(4.30) 


For the algorithm (4.30), we can consider the following basic objective functions. 


Maximization of the integral coefficient of synchronism 


K 




(4.31) 


2 = 1 


Maximization of the length of a synchronous domain 

= max >1, G TTIf} . (4.32) 

Note that, if necessary, one can expand the set of objective functions. Moreover, one can 


impose additional constraints on condition (4.30). 


A naive implementation of the algorithm (4.30) leads to difficulties of computational 


character. They are associated with combinatorial explosion. Let us illustrate this situation 
by an example of the full T-alphabet, IT"”^! = 17: 

N = 2, = 355 687 428 096 000, 

N = 3, = 126 513 546 505 547170 185 216 000 000, 

Af 3> 2, Curse of dimensionality! 


Currently, the author has developed a suboptimal version of algorithm (4.30). The 


main idea of the approach is as follows: the algorithm operates only with those symbols from 
the set V'" for which the elements of the matrix A"*" are close to the maximum value. Thus, 
this algorithm is free of the need of complete enumeration of the set however, it does 

not guarantee the global optimum of the function F^"^. Note that this topic is the subject 
of our current research. 
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RUB 



5. Sample 

Let us demonstrate the capabilities of the tools developed by an example of the analy¬ 
sis of hnancial time series. The object of analysis is the time series of exchange rates of some 
world currencies (US dollar [USD], Euro [EUR], Japanese Yen [JPH], Swiss Franc [CHE], 
and British Pound [GBP] against Russian ruble). The analyzed period is from 01.01.1999 
to 31.12.2014. 

Note that the analysis of the Generalized T-synchronization has applied value in the 
context of research in macroeconomics and stochastic hnancial mathematics. The original 
data are taken from the official web-site of the Gentral Bank of Russia (Bank of Russia, ex¬ 
change rates, www.cbr.ru/eng/). The length of the time series is iP = 3 985 samples. The ini¬ 
tial time series are shown in Figurej^ Note that these time series have also been studied from 
the viewpoint of complete T-synchronization (the analyzed period is from 01.01.1999 
to 31.03.2013) and TQ-complexity (2^ . 

We carried out a detailed analysis of generalized T-synchronization for the USD/EUR 
pair. For comparison, we also evaluated the characteristics of complete synchronization. The 
set was constructed for the objective function Fq"^. As a result, we obtained the following 
values of the integral coefficient of synchronism: 

YIG = 0.174492, Y|A = 0.219433, Y|G = 0.222948, (5.33) 


where (5®|G, (5®|A, and are complete, anti-, and generalized regimes of synchronization, 
respectively. 


The results (5.33) show that the regime of generalized T-synchronization is charac¬ 
terized by the maximum value of the integral coefficient of synchronism. At the same time, 
from the structural point of view, this mode is characterized by the combination of direct- 
and anti-synchronism regimes (see Figure]^. 


Taking account of the factor of generalized synchronization for the pair USD/EUR 
modihes the time structure of its T-synchronism (see Figure]^. 

Figure shows that taking account of the factor of generalized synchronization sig- 
nihcantly reduces the maximum length of desynchronous domains. In this case, the mean 
length of synchronous domains also decreases. 








Generalized synchronization of multidimensional chaotic systems... 


11 


TO T1 T2 T3N T3P T4N T4P T5N T5P T6S T6 T6L T7S T7 T7L T8N T8P 



Figure 6. Matrix A^; grey cells are T-symbols included in the set green and orange squares demonstrate 
the structure of the set in two cases of complete T-synchronization (see Figure]^; green and orange 
ellipses are key T-symbols in the set 



1 ^° \ 2 5 10 20 50 100 200 0 5 


10 


15 20 


25 


Figure 7. Spectral density of domains SD and SD; (a, b), (c, d), and (e, f) are complete, anti-, and generalized 
regimes of synchronization, respectively. 
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6. Conclusion 


We have proposed a new method to diagnose generalized synchronization in nonlinear 
multidimensional chaotic systems. The method allows one to explore and quantitatively 
evaluate the time structure of synchronization of chaotic oscillations in the so-called T- 
synchronization regime 16,17. The approach is based on the formalism of symbolic CTQ- 


analysis, which was proposed by the author in 18,19 


The method considered, which is based on the analysis of generalized T-synchronization, 
can be successfully applied to the study of multidimensional systems consisting of two or 
a greater number of coupled nonidentical oscillators, including multidimensional lattices of 
oscillators with arbitrary topology. The approach described can be applied to the analysis 
of experimental data, because it does not require any a priori knowledge of the system under 
study. 
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